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We study the problem of a Dirac field in the background of an Aharonov-Bohm flux 
string. We exclude the origin by imposing spectral boundary conditions at a finite 
radius then shrinked to zero. Thus, we obtain a behaviour of the eigenfunctions 
which is compatible with the self-adjointness of the radial Hamiltonian and the 
invariance under integer translations of the reduced flux. After confining the theory 
to a finite region, we check the consistency with the index theorem, and discuss 
the vacuum fermionic number and Casimir energy. 



1 Setting of the problem 

We study the Dirac equation for a massless particle in four dimensional Minkowski 
space, in the presence of a flux tube located at the origin, i.e., 

(i#-4)* = rl=VAi = -S(r)e z (1) 

r 

where n = S- is the reduced flux. 

n 

As the gauge potential is z— independent- equation (|l|) can be decoupled 
into two uncoupled two-component equations p by choosing: 
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In order to avoid singularities, we will consider that only r > r<j is accessi- 
ble, and take the limit ro — * 0, which is equivalent to having a punctured plane 
with the removed point corresponding to the string position. 

By taking A z = A r = 0, A$ = ^, for r > To, the Hamiltonian can be seen 
to be block-diagonal, with its two-by-two blocks given by 

/ ie* if> (d r ±B)\ i k 

H± -{-ze^ e (-d r ±B) ) B --r d °~r (3) 

It should be noticed that these two "polarizations" , which we will label 
with s = ±1, correspond to the two inequivalent choices for the gamma matri- 
ces in 2+1 dimensions. From now on, we will be working with s = 1 (the case 
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s = — 1 can be studied in a similar way, and explicit reference will be made to 
it whenever necessary). In this case, we can write: 

H + = ( °t n J . with L = -* e " i~ d r + B ) Li = * e ~ lf> ( d r + B) (4) 



(5) 



s L / 
and its eigenfunctions 
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Now, the two components in have different dependence. In order 
to make this fact explicit, and to discuss boundary conditions at r — r , we 
introducecl 

*E = ±[ e ~i! VlE ^\ L 1 = -d r + B L\=d r+ B (6) 
Vr \ e l 2 X iE (r,0) J 

so that 

L x ip XE = iE\iE l\xie = -iEipxE (7) 

We expand tpiE and X \e in terms of eigenfunctions of B, which are of the 
form: 

e n = e< n+ ?) e , with A» (r) = - + I - - , neZ (8) 

r 

once the condition has been imposed that ipE and xe in equation (|^) arc 
single- valued in 6. 
Thus, we have 

oc oo 

<piE(r,0)= £ / n (r)e<"+*) 9 XlJ3 (r,0) = £ Sn(r)e<"+^ (9) 
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which leads, for noninteger n (k = k + a, with k the integer part of k and a 
its fractionary part) to 

,t, < m W A„J„- K (|i?|r)+5 n J K - n (|£|r) 



^ (A, J n+ i_« (|£|r) - B n J K - n -i (\E\r)) 



e 



(10) 

(Of course, for integer k, a linear combination of Bessel and Neumann 
functions must be taken). Finally, for s = — 1, the upper and lower components 
of f £ interchange, and E — > —E. 
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2 Boundary conditions at the origin 



As is well known(a the radial Dirac Hamiltonian in the background of an 
Aharonov-Bohm gauge field requires a self-adjoint extension for the critical 
subspace n = k. In fact, imposing regularity of both components of the Dirac 
field at the origin is too strong a requirement, except for integer flux. Rather, 
one has to apply the theory of Von Neumann deficiency indices, which leads 
to a one parameter family of allowed boundary conditions, characterized by 

t lim (Mr)" +1 g n (r) sin + f ) = Hm (Mrp /„(r) cos ^ + |) (11) 

with v varying between —1 and (y — —a for s = 1; v = a — 1 for s = — 1). 
Here, parametrizes the admissible self-adjoint extensions, and M, a mass 
parameter, is introduced for dimensional reasons. Which of these boundary 
conditions to impose depends on the physical situation under study. 

One possibility is to take a finite flux tube, ask for continuity of both 
components of the Dirac field at finite radius and then let this radius go to zeroB 
Thus, one of the possible self-adjoint extensions is obtained, which corresponds 
to 9 = ^sgn(K). As pointed out, for instance, inp this kind of procedure 
leads to a boundary condition that breaks the invariance under k — ► k + n 
(neZ). Now, this is a large gauge symmetry, which of course is singular when 
considering the whole plane, but is not so when the origin is removed or, 
equivalently, the plane has the topology of a cylinder. 

To preserve the aforementioned symmetry we propose, instead, to exclude 
the origin, by imposing spectral boundary conditions of the Atiyah-Patodi- 
Singer(APS) typeP as defined inp at a finite radius tq and then letting tq — > 0. 

We consider the development in eq.(^) and, for s = 1, impose at r = tq : 

/„ (ro) = 0, for A„ (r ) < g n (r ) = 0, for A„ (r ) > (12) 

As is well known, imposing this kind of boundary condition is equivalent 
to removing the boundary, by attaching a semi-infinite tube at its position and 
then extending the Dirac equation by a constant extension of the gauge field, 
while asking that zero modes be square integrableO (except for A„ = 0, where 
a constant zero mode remains, with a nonzero lower component ). 

After using the dominant behaviour of Bessel functions for small argu- 
ments, and taking the zero radius limit, we have the following result for the 
eigenf unctions in eq. (|To|) : 

* f m V R ( J k+a -n{\E\r)e m9 
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n=k+l 



A„ 



Jn-k-a (\E\r) 



a > 



(13) 



Jn-fe-a (|^|r) e me \ 1 



i±fj n+1 ^ a (\E\r)e^+V e J "2 



a < - (14) 



Notice that our procedure leads precisely to a self adjoint extension satisfy- 
ing the condition of minimal irregularity (the radial functions diverge at r — > 
at most as r~ p , with p < ~). It corresponds to the values of the parameter : 

HI (15) 

As shown in^ 9 = ± ^ are the only two possible values of the parameter 
which correspond to having a Dirac delta magnetic field at the origin. More- 
over, this extension is compatible with periodicity in k. In fact, the dependence 
on k can be reduced to an overall phase factor in the eigenfunctions. 

As regards charge conjugation (^e —> cti^e*', k — > —k), it is respected by 
the eigenfunctions, except for a = 4. This is due to the already commented 
presence of a constant zero mode on the cylinder. However, for the repre- 
sentation s = -1 of 2 x 2 Dirac matrices, APS boundary conditions must be 
reversed, for A 0, as compared to (|l2|), since the operator B changes into 
— B. For A = the lower component will be taken to be zero at ro which, as we 
will show later, allows for charge conjugation to be a symmetry of the whole 
model. In this case, the resulting extension corresponds to 

e = { f ** a *\ ( i6) 

It is worth pointing that, for integer n — k, our procedure leads (both for 
s = ±1) to the requirement of regularity of both components at the origin. 



3 The theory in a bounded region 

From now on, we will confine the Dirac fields inside a bounded region, by 
introducing a boundary at r — R, and imposing there boundary conditions of 
the APS type, complementary to the ones considered at r = r$. For s = 1 

/„ (R) = 0, for A„ (JZ) > g n (R) = 0, for A„ (R) < (17) 
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We start by studying the zero modes of our theory, which are of the form 
* o (r,0)= e « "=-°° ^ J (18) 

It is easy to see that no zero mode remains after imposing the boundary 
conditions at both the internal and external radius, even without taking ro — > 
0. This is in agreement with the APS index theorem!] In fact, according to 
such theorem 

= A + b{r ) + b(R) (19) 

where rt + (n_) is the number of chirality positive (negative) zero energy solu- 
tions, A is the anomaly, or bulk contribution, and b are the surface contribu- 
tions coming from both boundaries Era 

b(R) = ~(h R - V (R)) b(r ) = ~(r)(r )-h ro ) (20) 

with rj(r) the spectral asymmetry of the boundary operator B and h r the 
dimension of its kernel. 

In our case, the boundary contributions cancel. As regards the volume 
part, it also vanishes for the gauge field configuration under study, and we 
have n+ — n_ = 0, which is consistent with the absence of zero modes. For 
s = — 1 both boundary contributions interchange, and identical conclusions 
hold regarding the index. 

The nonzero energy spectrum can be determined by imposing (for s = 1) 
the boundary conditions (|l7|) at r = R on the eigenfunctions in eqs. @ and 
@. Thus, one gets: 

j±h^L n >l I f±h=BA n >0 1 

E n i — { ,• , , , a > - E n i — { , H , , a < - (21) 

' \±^^-,n>-l ~2 ' |±J^i jn >o' 2 v ; 

where j Ut i is the i-th positive root of J v . The same spectrum results for s = — 1. 

For both s values, the energy spectrum is symmetric with respect to zero. 
This fact, together with the absence of zero modes results in a null vacuum 
expectation value for the fermionic charge □ 

(N) + = ~(n+-nJ)=0 (22) 

For the same reasons (iV)_ = 0, so the total fermionic number of the 
theory is null. 
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It is interesting to note that the contribution to the fermionic number 
coming from ro coincides, for each s value, with the result presented for the 
whole punctured plane ina(for details, setfcj), except for a = i. This last fact is 
associated with charge conjugation non invariance in each subspace. However, 
the sum of both contributions cancels for all a. 

We go now to the evaluation of the Casimir energy which, in the framework 
of the £ regularization, is given by 



where the parameter fi was introduced for dimensional reasons. Here, ( z ) is 
the so-called partial zeta function, defined as in referenceEH 
For the problem at hand, 



Notice that, as a consequence of the invariance properties of the imposed 
boundary conditions, the Casimir energy is both periodic in n and invariant 
under a — > 1 — a, as well as continuous at integer values of n. 

For any value of k, it is the sum of the energy corresponding to a scalar 
field in the presence of a flux string and subject to Dirichlet boundary con- 
ditions, plus a partial zeta coming, for fractionary k, from the presence of an 
eigenfunction which is singular at the origin or, for integer k, from the duplica- 
tion of Jo. Both contributions can be studied following the methods employed 
inllil So, here we won't go into the details of such calculation. 

The scalar field contribution presents a pole at z = —1, with an en- 
independent residue. However, the pole in the partial zeta appearing in (|23|),(p4[ 
has an a-dependent residue. Due to this fact, and the consequent need to in- 
troduce a-dependent counterterms, an absolute meaning cannot be assigned 
to the finite part of the Casimir energy. 

As a final remark, it is worth stressing that this approach to the problem 
of self-adjointness at the origin is, to our knowledge, the first proposal of a 
physical application of APS boundary conditions in this context. 

Contrary to the treatment of the origin, APS boundary conditions were 
imposed at the exterior boundary for merely formal reasons (the existence of 
an index theorem for this case). The vacuum energy under local (bag- like) 
external boundary conditions is at present under study. 





for oi>\ 
for a < \ 



(24) 
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